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Eulerian and Lagrangian tools are used to detect coherent structures in the velocity and 
magnetic fields of a mean-field dynamo, produced by direct numerical simulations of 
the three-dimensional compressible magnetohydrodynamic equations with an isotropic 
helical forcing and moderate Reynolds number. Two distinct stages of the dynamo are 
studied, the kinematic stage, where a seed magnetic field undergoes exponential growth, 
and the saturated regime. It is shown that the Lagrangian analysis detects structures 
with greater detail, besides providing information on the chaotic mixing properties of 
the flow and the magnetic fields. The traditional way of detecting Lagrangian coherent 
structures using finite-time Lyapunov exponents is compared with a recently developed 
method called function M. The latter is shown to produce clearer pictures which read- 
ily permit the identification of hyperbolic regions in the magnetic field, where chaotic 
transport/dispersion of magnetic field lines is highly enhanced. 
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1. Introduction 



The description of chaotic and turbulent flows by means of embedded coherent struc- 
tures is a topic of great interest in the study of transport and mixing in fluids, since these 
structures act as organizing units in the flow, defining attracting and repelling directions, 
transport barriers and regions of high or low dispersion of passive scalars. There is no 
standard way of defining what a coherent structure is, but from the Eulerian point of 
view, they are often defined based on some measure related to vorticity. An example is 



the highly popular Q-criterion, first introduced by Hunt et at. ( 1988 ) to identify vortex 



cores based on the difference between the rate of strain and vorticity. Some other criteria 



define coherent structures or vortices based on local pressure minima ( Jeong & Hussain 



1995 ) or on quantities involving the eigenvalues of the gradient tensor of the velocity 
field ( |Chong et al||1990[ |Zhou et a/.||1999[ |Chakraborty et aL||2005| |Varun et ai\20m\ . 
From a Lagrangian point of view, coherent structures are seen as material surfaces that 
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form the boundaries between regions of the flow with different behavior, such as vor- 
tex surfaces. They are found by following trajectories of fluid particles, while computing 



quantities such as the maximum rate of divergence of neighboring trajectories (Haller 
2001} |Shadden ef"aT|2005| ) or the arc-length of the trajectory ( |Madrid fc Mancho|2009f 



Lagrangian tools are naturally suited for unsteady flows, since they take into account the 
temporal variations of the vector field, not just instantaneous snapshots. For a recent list 
of applications, see |Peacock fc Dabiri (2010). 

Most works on Lagrangian coherent structures (LCSs) have focused on hydrodynamic 
turbulence, mainly in two-dimensions. A few papers have computed LCSs for three- 



dimensional magnetohydrodynamic (MHD) systems in the conservative (Leoncini et al. 
|2006[ ) and dissipative ( |Rempel et aT}|2011} |2012[ ) regimes. In the aforementioned dissi- 
pative cases, only velocity field (kinetic) structures were explored. Here, we expand our 
previous results by computing the kinetic and magnetic coherent structures in a MHD 
model of mean-field dynamo. Dynamo action consists in the amplification of magnetic 
field by the motion of an electrically conducting fluid, being the mechanism responsible 



for the equipartition-strength magnetic fields observed in planets and stars (Branden- 



burg fc Subramanian|2005[). Initially, a weak magnetic field B undergoes an exponential 



growth in the kinematic dynamo phase until B is strong enough to impact the fluid ve- 
locity u, and eventually the magnetic energy saturates. The saturation process is closely 
related to the suppression of Lagrangian chaos in the velocity field; a comparison be- 
tween the chaoticity of the velocity field during the growth and saturation phases of the 
dynamo has been performed in previous works ( Brandenburg et aL|1995 Cattaneo et al. 



1996 Zienicke et al. 1998). In this paper, the emphasis is on the detection of coherent 



structures and the transport of passive scalars and magnetic field lines in the transition 
from the kinematic to the saturated phase. Eulerian structures are detected using the 
Q-criterion and, for the detection of LCSs, the traditional technique of finite time Lya- 
punov exponents (FTLEs) is compared with the recently proposed function M ( Madrid 
& Mancho||2009|). 



Section [2] of this paper describes the dynamo model adopted. The numerical results 
are presented in section [3j where the Eulerian and Lagrangian coherent structures are 
computed for the velocity and magnetic fields. Some conclusions are given in section [4j 



2. The Model 

The model is the prototype of a 2 dynamo used by Brandenburg (2001 ), where a com- 
pressible isothermal gas is considered, with constant sound speed c s , constant dynamical 
viscosity p, constant magnetic diffusivity rj, and constant magnetic permeability po. The 
following set of compressible MHD equations is solved 



<%lnp + u-Vlnp + V-u = 0, (2.1) 
d t u + u-Vu = -Vp/p + J x B/p + {ji/p) (V 2 u + VV • u/3) + f , (2.2) 
d t A = u x B - np J, (2.3) 

where p is the density, u is the fluid velocity, A is the magnetic vector potential, J = 
V x B//io is the current density, p is the pressure, f is an external forcing, and Vp/p = 
c^Vlnp, where c 2 — jp/p is assumed to be constant. Nondimensional units are adopted 
by setting ki — c s = po = po = 1, where po = (p) is the spatial average of p and k\ 
is the smallest wavenumber in the box, which has sides L = 2ir and periodic boundary 
conditions. Thus, the time unit is (c s fci) _1 , space is measured in units of fc^ 1 , u in units 
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of c s , B in units of (popo) 1 / 2 ^, p in units of po and the unit of viscosity v = p/po 



and magnetic diffusivity 77 is c s jk\. Equations ( 2.1 )— ( 2.3 1 are solved with the PENCIL 
CODE [f] which employes an explicit sixth-order finite differences scheme in space and a 
third-order Runge-Kutta scheme for time integration. 

The initial conditions are Inp = u = 0, and A is a set of normally distributed, un- 
correlated random numbers with zero mean and standard deviation equal to 10~ 3 . The 
forcing function f is given by 

f (x, t) = Re{iVf k(t) exp[ik(t) ■ x + i(f)(t)]}, (2.4) 

where k(t) = (k Xl k y , k z ) is a time-dependent wavevector, x = (x,y,z) is position, and 
(f)(t) with \(f>\ < 7r is a random phase. On dimensional grounds the normalization factor is 
chosen to be N = /oC s (fcc s /<5t) 1 / 2 , where /q is a nondimensional factor, k — |k|, and St is 
the length of the integration timestep. We focus on the case where |k| is around kf = 5 
and randomly select, at each timestep, one of 350 possible vectors in 4.5 < |k| < 5.5. 
The operator fk is given by 



(2.5) 



ik x (k x e) — |k|(k x e) 
k ~ kV2(l- (k- e) 2 )/k 2 
where e is an arbitrary unit vector needed in order to generate a vector k x e that is per- 
pendicular to k. Note that |fk| 2 = 1 and the helicity density satisfies f • V x f = |k|f 2 > 0, 
which is an important condition for the production of a mean-field dynamo ( jMoffatt 



1978). The forcing function is delta-correlated in time, i. e., all points of f are correlated 



at any instant in time but are different at the next time step. Following Brandenburg 



(2001), the control parameters are set as fa = 0.07, v = -q = 0.002 and the numerical 
resolution is 128 3 . 



3. Results 

3.1. Mean-field dynamo 

1 /2 1/2 

Figure 1 shows the time series of B rms = (-B 2 ) (light line) and u rms = (u 2 ) (dark 
line), where B = |B| and u = |u|. During the first time units up to t ~ 150, the magnetic 
energy is too week to significantly impact the velocity field and u rms ~ 0.28, thus the 
Reynolds number is Re — u rms jvkj ~ 28. During this kinematic phase, B rms increases 
exponentially, with a growth rate 7 ~ 0.064 ± 2 x 10 -5 obtained from the fitted line 
(dashed line). After t ~ 150, u rms starts to decay due to the contribution of the Lorentz 
force (second term in the right side of Eq. (2.2 )). Eventually, the r.m.s. quantities saturate 
due to nonlinear effects, with u rms ~ 0.18 while the magnetic field reaches a super- 
equipartition value B rms ~ 0.37 > u rms . The arrows indicate the times t = 100 and 
t = 1700, respectively, which will be used later to represent the kinematic and saturated 
phases. In turnover time units (l/kfu rms ), the referred times are u rms kft ~ 140 and 
Urmskft ~ 1530, respectively, and the growth rate is j/urmskf ~ 0.046. 

During the kinematic stage, the magnetic field displays low-amplitude stochastic fluc- 
tuations, as shown in the upper panels of Fig. [2} As B rms grows, small-scale velocity and 
magnetic field fluctuations combine to produce a robust large-scale mean-field pattern 
(lower panels). The physics behind the rise of this mean-field is related to the so-called 
a-effect ( Moffatt||1978 ) and has been explored in this model by Brandenburg (2001). 



f http://pencil-code.googlecode.com 
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Figure 1 . Time series of B rms (light line) and u rms (dark line) of MHD dynamo simulations for 
n — v = 0.001. The arrows indicate the kinematic phase at t — 100 and the saturated nonlinear 
regime at t = 1700, respectively. The growth rate during the kinematic phase is 7 ~ 0.064. 



a) t = 100 








Figure 2. Intensity plot of magnetic field components at t = 100 (upper panel) and t = 1700 

(lower panel). 



3.2. Eulerian coherent structures 

Eulerian coherent structures can be extracted from the velocity field by decomposing the 
gradient tensor Vu as 



Vu = S + fJ, 



(3.1) 
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a) t=100 b) t=1700 




Figure 3. Eulerian coherent structures in the velocity field, detected by instantaneous 
isosurfaces of the Q-criterion. The isosurfaces are defined using 15% maximum Q. 

a) t=100 b) t=1700 




Figure 4. Eulerian coherent structures in the magnetic field, detected by isosurfaces of the 
Q-criterion. The isosurfaces are defined using 15% maximum Q. 



where S = ^ [Vu + ( Vu) ] and J~2 — ^ [Vu — (Vu) ] are the symmetric and antisym- 
metric parts of Vu, respectively. The symmetric part is the rate-of-strain tensor and 
the antisymmetric part is the vorticity tensor. One way to define an Eulerian coherent 
structure is by finding regions of u where vorticity dominates over strain, which can be 
measured by the Q-criterion (Hunt et al. 1988 Zhong et aZ.||1998 Haller|[2005 Lawson 
fc BarakospOlO ) 



(3.2) 



Thus, an Eulerian coherent structure or vortex is defined as a region where Q > 0. 

Figure [3] shows the isosurfaces of the Q _c riterion, us i n g 15% maximum Q (contour 
surfaces enclose high Q values). These plots are highly dependent on the threshold chosen 
for Q, but it is possible to sec that the fluid is more intermittent at the kinematic dynamo 
phase (t = 100) than after saturation (t = 1700), since in the right panel the coherent 
structures fill the space in a more homogeneous way. There are fewer regions for t = 100 
where Q is much higher than the average, thus the presence of fewer vortices for this 
threshold in Fig. |3]ja) than in[3]jb), where local values of Q are closer to the average Q. 
Figure [4] shows the corresponding plots of Q for the magnetic field, where the coherent 
structures represent magnetic vortices or current structures Brandenburg et al. ( 1996 1 . 
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Figure 5. Eulerian coherent structures in the velocity field, detected by the Q-criterion at 
t = 100 (left panel) and t = 1700 (right panel). 



In Fig. [5j intensity plots of the Q-criterion are shown for two-dimensional slices of the 
box at planes z — (upper panels) and x — (lower panels) at times t = 100 (left panels) 
and t = 1700 (right panels), respectively. Coherent structures with strong vorticity are 
observed as bright spots, such as the one highlighted by a box in Fig. [5ja). Notice that 
at t = 1700 a large number of bright spots is seen in the xy-plane, but they are rare 
in the yz-plane, revealing a preferential alignment of coherent structures in the vertical 
direction in the saturated regime. A similar plot is shown for the magnetic field in Fig.[6j 
where some of the same coherent structures found in the velocity field can be observed, 



reflecting the strong coupling between B and u in Eqs. (2.2) and (2.3) 



Although some coherent structures are clearly detected by this Eulerian technique, the 
Q-criterion relies on a user defined threshold to determine their boundaries. In order to 
precisely identify the boundaries and the main transport barriers in the flow, the next 
section proceeds with a Lagrangian analysis. 



3.3. Lagrangian coherent structures 

This section describes two tools that can be employed to define/detect coherent structures 
in the Lagrangian frame, the finite-time Lyapunov exponents and the recently proposed 
function M. 
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Figure 6. Eulerian coherent structures in the magnetic field, detected by the Q-criterion at 
t = 100 (left panel) and t = 1700 (right panel). 



3.3.1. Finite-time Lyapunov exponents 

In the Lagrangian point of view, coherent structures are seen as material surfaces 



around which trajectory patterns are formed. In Haller & Yuan (2000), these surfaces are 



simply called Lagrangian coherent structures and are distinguished from other material 
surfaces in that a LCS exhibits locally the strongest attraction, repulsion or shearing in 
the flow. Repelling LCSs are responsible for generating stretching, attracting LCSs for 
folding, and shear LCSs for swirling and jet-type tracer patterns ( |Haller|[2011 ). 

Attracting LCSs have commonly been associated with local maximizing curves (ridges) 
in the backward-time finite-time Lyapunov exponent (FTLE) field and repelling LCSs 
to ridges in the forward-time FTLE field ( Shadden et aZ.|2005j Green et aL|2007 Beron- 
Vera & Olascoaga 2010). There are limitations in such definition, as pointed out by Haller 
(2011) and Farazmand & Haller (20121, e. g., a ridge in the FTLE field may indicate 



the presence of a shear LCS or no LCS at all. But in general, ridges in the FTLE fields 
provide a good approximation to the true LCSs of the flow. 

Let D C R 3 be the domain of the fluid to be studied, let x(to) G D denote the position 
of a passive particle at time to and let u(x, t) be the velocity field defined on D. The 
motion of the particle is given by the solution of the initial value problem 



— =u(x(t),t), x(t )=x . 



(3.3) 
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Let the flow map for u be defined as (f>l°. +T : x(io) x(4o+t). The deformation gradient 
is given by J = d</>j° +T (x)/dx and the finite-time right Cauchy-Green deformation tensor 
is given by A = J J. Let Ai > A 2 > A 3 be the eigenvalues of A. Then, the finite-time 
Lyapunov exponents or direct Lyapunov exponents of the trajectory of the particle are 
defined as 



a*° +T (x) = ~ln^/\i, i = 1,2,3. (3.4) 

m 

A positive 0\ is the signature of chaotic streamlines in the velocity field, being a 



measure of the stretching of fluid elements (although it also incorporates shear ( Haller 



2011)). In this work, the deformation gradient is computed with second order centered 



finite-differences. 
3.3.2. Function M 



Madrid & Mancho (2009) proposed a function to define "distinguished trajectories" 
(DTs), which are a generalization of the concept of fixed points for aperiodically time- 
dependent flows. In stationary flows, hyperbolic fixed points are responsible for particle 
dispersion and nonhyperbolic fixed points for particle confinement. Invariant stable and 
unstable manifolds of hyperbolic fixed points are barriers to transport and divide the 
phase space in regions with qualitatively different behaviors. The proposed function, 
named function M, can reveal both hyperbolic and nonhyperbolic flow regions of time- 
dependent flows. Moreover, M is also useful to detect the stable and unstable manifolds 
of distinguished hyperbolic trajectories (DHTs), defined as the set of points such that 
trajectories passing through these points at t = to will approach the DHTs at an expo- 
nential rate as time goes to infinity or minus infinity, respectively ( jBranicki et «Z.|2011 ). 



The stable and unstable manifolds of DHTs correspond to the repelling and attracting 
Lagrangian coherent structures, respectively, as defined in the section |3.3.1| 



Consider the system given by Eq. (3.3 ), where x = {x±, X2, x^}. For all initial conditions 



x in D at a given time to, let us define the function M(x. ,t ) : (D,t) — > K as 



rto+r ( 3 \ 1 I 2 

M(xo,t )r= [y2(dXi(t)/dt) 2 ] dt. (3.5) 

Jt -T \. = 1 J 

Thus, the function M is a measure of the arc length of the curve traced by Xo- Local 
minima of M represent trajectories that "move less" , being related either to hyperbolic 
or nonhyperbolic DTs. The manifolds of DHTs are also visible in the M field, since one 
expects a sharp distinction in the lengths of trajectory curves for particles in regions with 
different behaviors, separated by stable and unstable manifolds, as noted by |Mendoza fc| 



Mancho (2010). The technique has been successfully applied to the detection of DTs and 
manifolds in oceanic (|Mendoza fc Mancho||2010 Mendoza et a^.pOlO ) and stratospheric 



(de la Camara et al. 2Q12J flows. 



3.3.3. Velocity field structures and chaotic mixing 

The FTLEs are computed from a series of fully 3D snapshots of the velocity field 
taken at different times from to to to + t. Linear interpolation in time and third-order 
splines in space are used to obtain the continuous vector fields necessary to obtain the 
particle trajectories. Figure [7] depicts the probability distribution functions (PDFs) of 
the three FTLEs at t = 100 (left) and t = 1700 (right) computed for 64 3 particle 
trajectories from Eq. ( |3.4[ ) with a value of r corresponding to 9 turnover time units, 
where u rms ~ 0.28 for the kinematic phase and u rms ~ 0.18 for the saturated regime. 
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Figure 7. Probability density functions (PDFs) of the FTLEs of the velocity field at the kine- 
matic (t — 100, left panel) and saturated (t = 1700, right panel) regimes. The solid line represents 
cti, the dotted line, 02, and the dashed line, 0-3. 



Therefore, r = 9/(kfU rms ) ~ 6.4 time units for the kinematic phase and r ~ 10 time 
units for the saturated phase. One can see a clear reduction of Lagrangian chaos in 
the velocity field at t — 1700, with the PDF of a x being shrunk and shifted to the 
left. There are also fewer regions with two or three positive exponents. Overall, chaotic 
mixing is diminished due to the growth of B rms and the action of the Lorentz force. The 
asymmetry in the distributions is typical of heterogeneous mixing, where both regular 
and irregular trajectories coexist (in finite-time), which means that trajectories cannot 
uniformly sample the phase space (see, e. g., |Beron-Vera fc Olascoaga (2010)). 

From Fig. [7] it is clear that most trajectories display two positive Lyapunov exponents. 
Zel'dovich et al. ( 1984 ) and Chertkov et al. ( 1999 1 state that in such a case, the total 
magnetic energy in a kinematic dynamo should behave as B 2 oc exp[(cri — cr 2 )^, therefore, 



one has for the growth rate 



7 = 



d\n(B 2 ) 1/2 

dt " 
0.339 and (ct 2 ) 



dln{exp[(cri - a 2 )t/2}} _ a x - ct 2 
dt ~ 2 ' 



(3.6) 



At t = 100, (cti) - 0.339 and (ct 2 ) ~ 0.143, which from Eq. (|3_6|) provides 7 = 0.098 (or 
7 = 0.07 in dimensional units), which agrees to within an order of magnitude with the 
fitted value 7 ~ 0.064, given in Fig. [I] 

The remainder of this paper focuses on the backward-time maximum FTLE field, 
since they reveal the attracting LCSs, which correspond to structures seen using flow 
visualization in experiments (Green et al. 2007). Figure [8] shows the backward-time 
maximum FTLE field computed for r corresponding to 9 turnover time units at to = 100 
(left) and to = 1700 (right) from a grid of initial conditions with 512 x 512 particles. The 
bright lines represent the attracting LCSs. While the LCSs at to = 100 reveal no preferred 
direction, consistent with an isotropic forcing, at to = 1700 there is a clear vertical 
alignment of LCSs in the iry-slice (Fig. [8](b)) . This is due to the super-equipartition 
magnetic field at to = 1700, which develops a large-scale vertical pattern in this plane 
(see Fig. [2]), affecting the alignment of velocity field vectors. 

A comparison between Figs. [5] and [8] shows that the FTLE field provides a clearer 
depiction of coherent structures, with finer details and more precise detection of structure 
boundaries. Moreover, some coherent structures are only apparent in the FTLE field, such 
as the large eddy indicated in Fig. ^d) . 

From our experience, one of the problems with FTLE plots in turbulent flows is that 
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Figure 8. Attracting Lagrangian coherent structures in the velocity field, given by the 
backward-time FTLE at to = 100 (left panel) and to — 1700 (right panel). 



pictures usually become increasingly complex for larger r, with material lines "growing" 
and filling the entire phase space. In that sense, it is easier to use function M to detect 
the main coherent structures of the flow. Figure [9] is a plot of function M with r = 9 
in turnover time units, equivalent to Fig. [8j For the kinematic phase at to = 100 (left 
panel) the eddies are clearly identified as closed regions with distinguished shades of gray. 
At to = 1700 (right panel) the borders between regions are not so sharp and there are 
wide smooth regions in the flow. Smoothness in the M field indicates that trajectories in 
those regions do not reach nearby hyperbolic regions during ±r turnover time units, since 
hyperbolic trajectories are the ones responsible for dispersion and for producing sharp 
changes in M (Mendoza & Mancho 2010). For larger r, the boundaries become sharper 
and more foldings of manifolds are seen, but we keep r = 9 in all our pictures to facilitate 
the comparison between both methods in different regimes. Overall, the function M seems 
to be less sensitive to the choice of r than the FTLE. Figure [9] corroborates with Fig. [7] 
in revealing that there is less chaotic mixing after the nonlinear saturation of B rms . 



3.3.4. Magnetic field structures and transport of field lines 

Our simulations reveal that the magnetic field displays smooth and complex regions 
(see Fig. |6| . If one applies the Lagrangian techniques discussed in the previous section to 
the magnetic field, the identification of magnetic LCSs provides the main barriers to the 



Coherent structures in a nonlinear dynamo 



11 
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Figure 9. Lagrangian coherent structures in the velocity field, given by the function M at 
to = 100 (left panel) and to = 1700 (right panel). 



transport of field lines, a topic of great interest in magnetic reconnection studies ( |Evans 
et a/.||2004| |Grasso et aL|2010| |Borgogno eFa^On! |Yeates fc Hornig||2011[ ) . 



To obtain the magnetic LCSs, the magnetic field at a fixed dynamic time to is used 
and the maximum FTLE field is computed by integrating 



f/x 

ds 



B(x(s),< ), x(s )=x , 



(3.7) 



where the parameter (position) s along the field line is seen as an effective time, or field- 



line-time (Borgogno et al. 20111. The flow map for B is defined as 



s 



x(so) 



x(s + t). Equation ( |3.7[ ) is integrated from s to s + r with t fixed at to- Lagrangian 
chaos in the magnetic field is responsible for the transport of magnetic field lines between 
different regions of the box. Here, the term "transport" is used to refer to motion of field 
lines in field-line-time, not in dynamic time. Therefore, the maximum FTLE provides a 
measure of the exponential separation between two neighboring field lines after a finite 
field line time r, i. c.. after a finite distance along the field line. 



Figure 10 shows the backward-time maximum FTLE field for the kinematic (left panel) 



and saturated (right panel) regimes. The high-intensity lines represent attracting mag- 
netic LCSs which act as barriers to field line transport. No transport of magnetic field 
lines occurs across invariant LCSs and large-scale transport is possible only through 
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t = 100 t = 1700 




y 



Figure 10. Attracting Lagrangian coherent structures in the magnetic field, given by the 
backward-time maximum FTLE at t = 100 (left panel) and t — 1700 (right panel). 



homoclinic and heteroclinic crossings of attracting and repelling LCSs, where a lobe dy- 



namics mechanism takes place ( Grasso et aL|2010 


Borgogno et al. 2011 


Yeates & Hornig 


2011| Rempel et al. 


2012). Both FTLE fields are obtained by fixing t 


ic evolution (dy- 



namic) time (to 



t = 9/B rms , where B rms 



100 for the left panel and to = 1700 for the right panel) and setting 
0.014 for t = 100 and B rms = 0.37 for t = 1700. In the 
kinematic regime (to — 100) the LCSs display no preferred direction, and randomly fill 
the simulation box. Note that, at least for this value of r, it is difficult to identify the 
coherent structure marked in the box in Fig. [To{a) due to the many foldings of attracting 
lines. After growth and saturation of B rms (to = 1700), the randomness of field line 
orientation is diminished and there is a preferential direction of alignment of field lines 
which, as mentioned before, directly affects the velocity field. The growth of B rms is also 
reflected in the presence of many smooth regions in Figs. [Toj^c) and (d). Consequently, 
there is less chaotic mixing of field lines. 

Once again, to obtain a clearer picture of magnetic coherent structures, we plot in Fig. 

9/B rms . It is easier to spot coherent structures from this field. 
- a). 



11 



the function M for r 



such as the one in the box in Fig. 11 a). Function M seems to be better than the FTLE 
field in highlighting the main transport barriers, filtering out spurious lines that are not 
so important for mixing ( Mendoza fc Mancho||2010 ). 

As mentioned before, another feature of function M plots is that they provide both the 
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stable and unstable manifolds of DHTs in the same picture. In order to illustrate this 
feature, three distinct regions are marked in Fig. [TT^d) . Regions A and B are located 
in smooth parts of the M field and region C in a region where manifolds are crossing. 
Smoothness of M in regions A and B indicates that initial conditions in these regions do 
not perceive nearby hyperbolic regions for f G (to - T i to + T ) ( Mendoza &; Mancho|2010 ). 
An enlargement of region C is shown in Fig. [l2j where the presence of manifolds indicates 
that field lines in this region either were dispersed in to — r or will disperse in to + r. We 
define three sets of initial conditions inside the small white squares A, B and C in Fig. 
|Tl}d ), with each square containing 25 initial conditions. The result of integrating Eqs. 
(3.7) with each set of initial conditions for r = 9/B rms field-line-time units is shown 
in Fig. 13 Figures [13^ a) and (b) show the trajectories of initial conditions in regions A 
and B, respectively, where it can be seen that all magnetic field lines stay close to each 
other, forming a magnetic flux tube that is not dispersed in this field-line-time interval. 
The apparent discontinuities in field lines are due to the periodic boundary conditions. 
In Fig. |l3"{c) the trajectories of initial conditions at region C are shown and one can see 
that there is great chaotic dispersion of field lines due to the crossings of manifolds in 
this region. 

We conclude that function M can efficiently detect transport barriers and dispersion 
regions in a magnetic field. 
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Figure 13. Magnetic field lines produced by advecting a small blob of initial conditions in the 
magnetic field at t = 1700. (a) The initial blob is located at A, the point inside a magnetic 
vortex in Fig. JTTJd) ; (b) the initial blob is located at B, the point in a smooth region of Fig. 
|ll[ d); (c) the initial blob is located at C, the point at a crossing of manifolds in Fig. |ll[ d). 
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4. Conclusions 

Magnetohydrodynamic coherent structures have been identified in direct numerical 
simulations of a nonlinear dynamo. It was shown that both Eulerian and Lagrangian tools 
are able to extract vortices from velocity and magnetic field data. Although the Eulerian 
tool adopted is less computationally expensive, Lagrangian plots show finer details and 
can better locate the boundaries of vortices. In addition, the Lagrangian analysis provides 
important information about the mixing properties of the flow. Regarding the numerical 
tools employed to detect Lagrangian coherent structures (LCSs) , the function M seems 
to be less sensitive to the choice of the integration time r in comparison to the maximum 
finite-time Lyapunov exponent (FTLE). Thus, pictures obtained with the FTLE can 
become increasingly "noisy" with increasing r due to the complex folding of material 
lines. Although function M provides "cleaner" pictures, the manifolds (transport barriers) 
are often not as clearly traced as in a FTLE field. Both tools reveal the strong impact of 
the magnetic field on the mixing properties of the velocity field when the system moves 
from the kinematic to the saturated dynamo phases. After the appearance of a strong 
mean-field, the kinetic and magnetic coherent structures are shown to align in a preferred 
direction, revealing the anisotropy developed in the vector fields. 

Function M is also shown to be useful to detect manifolds of hyperbolic trajectories in 
the magnetic field, where intense transport of magnetic field lines takes place, a feature 
that can be further explored to study magnetic reconnection phenomena in plasmas. In 
relation to this, Lagrangian coherent structures in photospheric velocity fields have been 



shown to be associated with quasi-separatrix layers in the magnetic field (Yeates et al. 



2012), which are regions of strong gradients in stretching and squashing of magnetic flux 



tubes, being identified as the preferential regions for magnetic reconnection (Dcmoulin 



2006 Santos et al. 2008). Magnetic reconnection is an important phenomenon in non- 



linear dynamos, since it is believed that it can reduce the backreaction of the Lorentz 
force on the velocity field ( Blackman|1996 ). Essentially, turbulent motions can cause the 
stretching, twisting and folding of weak magnetic field lines in such a way as to pro- 
duce the growth of magnetic flux. After the magnetic field reaches equipartition with 
the velocity field, the field lines can restrict fluid motions and transport of material is 
significantly reduced. This suppression of motions may also inhibit the dynamo. How- 
ever, if there is rapid reconnection between magnetic flux tubes, this could prevent the 
tube from backreacting. For other works on the role of magnetic reconnection in dynamo 
models, see Archontis et al. (2003) and |Baggaley et al. (2009). 
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